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We study the space-time evolution of the fine structure constant, a, inside evolving spherical 
overdensities in a lambda-CDM Friedmann universe using the spherical infall model. We show that 
its value inside virialised regions will be significantly larger than in the low-density background 
universe. The consideration of the inhomogeneous evolution of the universe is therefore essential for 
a correct comparison of extragalactic and solar system limits on, and observations of, possible time 
variation in a and other constants. Time variation in a in the cosmological background can give 
rise to no locally observable variations inside virialised overdensities like the one in which we live, 
explaining the discrepancy between astrophysical and geochemical observations. 
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Recent observations of small variations in relativistic 
atomic structure in quasar absorption spectra Q sug- 
gest that the fine structure constant a, was smaller at 
redshifts z = 1 — 3 than the current terrestrial value 
a = 7.29735308 x 10" 3 , with Aa/a = {a(z)-a }/a = 
—0.543 ± 0.116 x 10~ 5 . Several theories of varyinga 
have been proposed to investigate the implications [2|. 
Time-varying a requires a scalar field that couples to 
electromagnetically-charged matter. Variations in a, al- 
lowed by the conservation of energy and momentum will 
then depend on the expansion of the universe and the 
evolution of the electromagnetically coupled matter. The 
inclusion of electroweak or grand unification will create 
more complicated consequences and constraints . Vari- 
ations in a due to perturbation in the matter fields were 
first studied in ||, using the linear theory of cosmologi- 
cal perturbations. It was shown there, that perturbations 
in a will grow during the matter-dominated epoch and 
will be constant or decay during the other eras. In Q 
the effects of a on gravity were also analysed. However, 
any effects on the time and spatial evolution of a due to 
cluster and galaxy formation has been ignored in the lit- 
erature. This is a major weakness and, as a result, past 
attempts to confront observations on extragalactic scales 
with lab or solar system constraints on a variation are 
all questionable. This letter describes the first attempts 
to follow the inhomogeneous evolution of a during the 
development of non-linear cosmic structures. We com- 
pare the evolution of a in the background universe and 
in overdense regions of the universe. 

We will study variations in a in the Bekenstein- 
Sandvik-Barrow-Magueijo (BSBM) theory 0, which as- 
sumes that the total action of the Universe is given by: 
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In the BSBM varying-a theory, the quantities c and 
h are constants, while e varies as a function of a real 
scalar field tp, with e = e e^ . — ^d^tpd^ip, uj 

is a coupling constant, and £ em = —\f l ivf iiV - The 
gravitational Lagrangian is the usual C g = — 



with R the Ricci curvature scalar. C m is the matter 
fields Lagrangian. Defining an auxiliary gauge poten- 
tial by a M = eAfi and a new Maxwell field tensor by 
ffj,u = e-F/ii/ = d^o-v — d„a^, the covariant derivative takes 
the usual form, = + ieoa fl . The fine structure con- 
stant is then a = a^e 2 ^ with ao the present value here. 

The background universe will be described by a flat, 
homogeneous and isotropic Friedmann metric with ex- 
pansion scale factor a(t). Varying the total Lagrangian 
we obtain the Friedmann equation (h = c = 1) for a 



universe containing dust, of density p„ 
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and a cos- 



mological constant A with energy-density p A = A/ (8t:G): 
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where H = a/ a is the Hubble rate, 



, and 



A/- — 2 ' 

C = £em/p m is the fraction of the matter which car- 
ries electric or magnetic charges. The value (and sign) 
of C will depend on the nature of dark matter: £ w 10~ 4 
for neutrons or protons but £ = — 1 for superconducting 
cosmic strings |4|. The scalar-field evolution equation is 
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Variations in a, due to the formation of overdense 
spherical regions, can be studied using the spherical infall 
model [lfH. We model the spherical overdense inhomo- 
geneity by a closed Friedmann metric and the background 
universe by a spatially flat Friedmann model. Each has 
their own expansion scale factor and time, which are 
linked by the condition of hydrostatic support at the 
boundary. This approach is equivalent to study the effect 
of perturbations to the Friedmann metric by considering 
spherically symmetric regions of different spatial curva- 
ture in accord with BirkhofF's theorem. The density per- 
turbations need not to be uniform within the sphere: any 
spherically symmetric perturbation will evolve within a 
given radius in the same way as a uniform sphere con- 
taining the same amount of mass. Clearly, this model 
ignores any anisotropic effects of gravitational instability 
or collapse. Similar results could be obtained by perform- 
ing the analysis of the BSBM theory using a spherically 
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symmetric Tolman-Bondi metric for the background uni- 
verse with account taken for the existence of the pressure 
contributed by A and the ip field. In what follows, there- 
fore, density refers to mean density inside a given sphere. 

Consider a spherical perturbation with constant inter- 
nal density which, at an initial time, has an amplitude 
Si > and \Si\ *C 1. At early times the sphere expands 
along with the background. For a sufficiently large Si, 
gravity prevents the sphere from expanding. Three char- 
acteristic phases can then be identified. Turnaround: 
the sphere breaks away from the general expansion and 
reaches a maximum radius. Collapse: if only gravity is 
significant, the sphere will then collapse towards a cen- 
tral singularity where the densities of the matter fields 
would formally go to infinity. In practice, pressure and 
dissipative physics intervene well before this singularity 
is reached and convert the kinetic energy of collapse into 
random motions. Virialisation: dynamical equilibrium 
is reached and the system becomes stationary: there are 
no more time variations in the radius of the system, R, 
or in its energy components. In the spherical collapse 
model, due to its symmetry, the only independent coor- 
dinates are the radius of the overdensity and time. Also, 
as is standard practice when using this model, we con- 
sider that there are no shell-crossing which implies mass 
conservation inside the overdensity and independence of 
the radius coordinate |lOj| . The equations can then be 
written ignoring the spatial dependence of the fields (but 
still including an equation for the evolution of the radius) . 
Hence, the evolution of a spherical overdense patch of 
scale radius R(t) is given by the Friedmann acceleration 
equation: 

3R = -iirGR (p cdm (I + |C| e" 2 ^) + 4p^ - 2p A ) (4) 

where p cdm is the total density of cold dark matter and 
baryons in the cluster, ip c is the homogeneous field in- 
side the cluster and we have used the equations of state 

P*P C = Pip c > P cdm = and P A = "^A- In tne duster; tne 
evolution of ip c and p cdm is given by 

</>' c + 3 (RIB) i> c = -2e- 2 HpcdJ" (5) 

and p c( j[ m cx R~ 3 . The cluster will virialise when 

T mr = -0.5f/ G + U A - 2U^ c (6) 

where Uq = — ^GM 2 R~ 1 is the potential energy as- 
sociated with the uniform spherical overdensity, U A = 
— ^irGp A MR 2 is the potential associated with A, and 
U$ = —^GMAI^ c R~ A is the potential associated with 
ip c . T V i r is the kinetic energy, and M = M c & m + M$ is 
the cluster mass, with M cdm = ^-p cdm (l + ICIe -2 ^)^ 3 
and Mjj! c — ^j-p^ R 6 ', Using the virial theorem © and 
energy conservation when the cluster virialises gives 

0.5U G + 2U A - \ z = Zv = U G + U A + U.^ \ z = zta (7) 



where z v is the redshift of virialisation and z ta is the red- 
shift at the turn-around of the over-density at its maxi- 
mum radius, when R = R max and R = 0. 

The analysis in this letter (which reduces completely 
to the linear perturbation analysis in the small time limit 
when all inhomogeneities are small 0) is valid for spher- 
ically symmetric perturbations right into the non-linear 
regime, including turnaround and collapse. Spatial gra- 
dients in pressure (and hence in the scalar field) are neg- 
ligible on large scales, exceeding the Jeans length. The 
scalar field is indeed slightly inhomogeneous but the in- 
clusion of a varying a at a level consistent with obser- 
vation 0, (a/ao ~ 10~ 6 H), does not affect the overall 
expansion of the universe or the dynamical collapse of the 
overdense regions. The energy density associated with ip 
is always a negligible contribution to the energy content 
of the universe, both in the cluster and in the background, 
if we are far from the initial or collapse singularities p|. 
After the fluctuation virialises and attains gravitational 
equilibrium this assumption may no longer be true, but 
our analysis like that of all other exact studies of galaxy 
formation breaks down at that stage because the result- 
ing disk will be controlled by pressure and gradients. 

The behaviour of the fine structure constant during 
the evolution of a cluster can now be obtained by numer- 
ically evolving the background eqs. (|2I>- ((31 and the clus- 
ter eqs.Q-JHJ until the virialisation condition @ holds. 
Since the Earth is in a virialised overdense region, the 
initial conditions for tp are chosen so as to obtain our 
measured lab value of a at virialisation a c (z v ) — ao and 
to match the latest observations 1} for overdense regions 
at 3 > \z — z v \ > 1. Since z v is uncertain, we have cho- 
sen representative examples with virialisation over the 
range < z v < 10. This is shown in Fig. |3 where 
the clusters have different initial conditions in order to 
satisfy the constraints described. This is just an exam- 
ple, since in reality, the initial condition for tp needs to 
be fixed only once, for our Galaxy. Hence, a in other 
clusters will have a lower or higher value (with respect 
to ao) depending on their 

FigH After viriali- 
sation occurs the cluster radius becomes constant; time 
and space variations in a are suppressed, and a becomes 
constant. If there were any variations in a after viriali- 
sation, the energy and radius of the cluster would need 
to vary in order to conserve energy and this is inconsis- 
tent with virialisation. This phenomenon is not included 
in Fig. ^ since we did follow the evolution to virialisa- 
tion with a many-body simulation which would need to 
include the fluid equations that describe the pressure in- 
side the cluster. In our simulation, the virial condition 
is a 'stop condition' and so we just observe the typical 
behaviour of the cluster's collapse as R — > 0. Clearly, 
collapse will never occur in practice; dissipative physics 
will eventually intervene and convert the kinetic energy of 
collapse into random motions. In addition to the station- 
arity condition that occurs when the cluster virialises, it 
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can be seen from FigJ^that in all cases the variation in 
a since the beginning of the cluster formation is of order 
10~ 5 , and numerical simulations give a/a ~ 10 _22 s _1 . 
If variations in a are so small for such a wide range of 
virialisation redshifts, we can assume that the difference 
between the value of a c at z v and at z = will be neg- 
ligible. Therefore it is a good approximation to assume 
that the time-evolution of both a c and of the cluster will 
cease after virialisation. Although this is not necessarily 
true (the cluster could keep accreting mass), it is a good 
approximation in respect of the evolution of a, especially 
for clusters which have virialised at lower redshifts. 
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FIG. 1: Evolution of a in the background (dashed line) and 
inside clusters (solid lines) as a function o/log(l + z). Ini- 
tial conditions were set to match observations of a variation 
in ref. jj] • Two clusters virialise at different redshifts, one 
of them in order to have a c (z v = 1) = Qo. Vertical lines 
represent the moment of turn-around. 



x 10~ 6 




FIG. 2: Evolution of Act/a in the background (dashed lines) 
and inside clusters (solid lines) as a function o/log(l + z). 
Initial conditions were set to match observations of a vari- 
ation in ref. Four clusters that virialised at different 
redshifts. All clusters were started so as to have a c (z v ) = Go- 



lf we could measure a inside virialised overdensities 
and the corresponding value in the background, at the 
time of their virialisation, then Fig [21 and 01 would give 
us the evolution of a v or its time shift, Aa v /a v , as a 
function of z v . Those figures show us that differences in 
between the background and the overdensities increase 
as z v — > 0. This is due to the earlier freezing of the 
value of a at virialisation, and to our assumption that 
we live in a AC DM universe. At lower redshifts, spe- 
cially after A starts to dominate, variations in a in the 
background are turned off by the accelerated expansion. 
However, the value of a in the collapsing cluster will keep 
growing until virialisation occurs. Numerical simulations 
give Aa/a(z = 2) = —6 x 10~ 6 for the background and 
Aa/a(z v = 0.13) = —1.5 x 1CP 7 in a cluster. At higher 
redshifts, z>l, both «{, and a c evolve in expanding en- 
vironments: their increase is logarithmic in time before 
A starts to dominate, so the difference between them will 
be much smaller. 
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FIG. 3: Plot of a as a function of\og(l + z v ), at virialisation. 
Clusters (solid line), background (dashed line). 

Spatial variations in a can be tracked using a 'spa- 
tial' density contrast, defined by Sa/a = [a c — ab]/a&. 
A plot of the spatial inhomogeneities in a with respect 
to the matter density inhomogeneity, Sp/p, at virialisa- 
tion is shown in Fig. Note that 6a/ a grows in pro- 
portional to the density contrast of the matter inhomo- 
geneities (tx i 2 / 3 when both are small during dust dom- 
ination ^5])- In a AC DM model, the density contrast, 
A c = p cdm (z v ) / p b (z v ) increases as the redshift decreases. 
For high redshifts, the density contrast at virialisation be- 
comes the asymptotically constant in standard (A = 0) 
CDM, A c « 178 at collapse or A c ps 148 at virialisation. 
This is another reason why at lower redshifts the differ- 
ence between a v and at, increases. Trends of variation 
in a can be then predicted from the value of the matter 
density contrast of the regions observed. Useful fitting 
formulas for the dependence of a variation on Sp/p and 
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FIG. 4: Plot of Aa/a as a function o/log(l + z v ), at viriali- 
sation. Clusters (solid line), background (dashed line). 

the scale factor a = (1 + are: 

Sa/a = (5.56 - 0.7(7* +0.078(7 + 0.00352(7*) x 10~ 6 
Sa/a = (5.37 4- 0.373(9* — 0.276** + 0.00761*) x 1(T 6 
Aa/a = -(2.47- 1.81cr* +0.59(7 - 0.094(7*) x 10~ 6 

where 9 = — /A Cij , a = a/a v . A Cv and a v are 'input' 
parameters defined when a{z v ) = a$ . 
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FIG. 5: Variation of Sa/a with Sp/p at the cluster viriali- 
sation redshift. The evolution of a inside the clusters was 
normalised to satisfy the latest time-variation observational 
results, and to have a c (z v ) = ao for z v — 1. 

The evolution of a = ab in the background and in- 
side clusters depends mainly on the dominant equation 
of state of the universe and the sign of the coupling con- 
stant — , which is determined by the theory and the dark 
matter identity. Here, we have assumed that lu = 1, so all 
the dependence is in £. As shown in refs. Q, ab will be 
nearly constant for accelerated expansion and also dur- 
ing radiation-domination far from the initial singularity 



(where the kinetic term, p^, can dominate). Evolution 
of a will occur during the dust-dominated epoch, where 
a increases logarithmically in time for £ < 0. When £ is 
negative, a will be a growing function of time but will 
fall for £ positive p|. Inside clusters a c will have the 
same dependence on £ as it has in the background. The 
sign of C is determined by the physical character of the 
dark matter that carries electromagnetic coupling: if it 
is dominated by magnetic energy then £ < 0, if not then 
C > 0. Our numerical simulations were performed assum- 
ing C = —2 x 10~ 4 . If we had chosen £ to be positive we 
would find that ab would decrease as steeply as z — > 0. 
We choose the sign of £ to be negative so a is a slowly 
growing function in time during dust domination. This is 
done in order to match the latest observations which sug- 
gest that a had a smaller value in the past Q]. However, 
we know that on small enough scales the dark matter will 
become dominated by a baryonic contribution for which 
C > 0. This will create distinctive behaviour and will be 
investigated separately. Generally, past studies of spa- 
tially homogeneous cosmologies have effectively matched 
the value of ab with the terrestrial value of a measured 
today. However, it is clear that the value of the fine struc- 
ture constant on Earth, and most probably in our local 
cluster, will differ from that in the background universe. 
This feature has been ignored when comparing observa- 
tions of a variations from quasar absorption spectra 
with local bounds derived by direct measurement or 
from Oklo and long-lived /?-decayers @ . A similar un- 
warranted assumption is generally made when using solar 
system tests of general relativity to bound possible time 
variations in G in Brans-Dicke theory |ll| : there is no 
reason why G/G should be the same in the solar system 
and on cosmological scales. Since any varying-constant's 
model require the existence of a scalar field coupled to the 
matter fields, our considerations apply to all other theo- 
ries besides BSBM and to variations of other 'constants' 
[l2|. Note that this feature is much less important when 
using early universe constraints like the CMB or BBN 
|l3l | , since small perturbations in a will decay or become 
constant in the radiation era [5j . 

In summary, using the BSBM theory we have shown 
that spatial variations in a will inevitably occur because 
of the development of large density inhomogeneities in 
the universe. This was first shown in the linear regime, 
when perturbations are small , and then £s tracks 
during the dust-dominated era on scales smaller than 
the Hubble radius. We have used the spherical collapse 
model to study the space-time variations in a in the non- 
linear regime. Strong differences arise between the value 
of a inside a cluster and its value in the background 
and also between clusters. Variations in a depend on 
the matter density contrast of the cluster region and the 
redshift at which it virialised. If the overdense regions 
are still contracting and have not yet virialised, then the 
value of a within them will continue to change. Varia- 
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tions in a will cease when the cluster virialises so long 
as it does so at moderate redshift. This leads to larger 
values of a in the overdense regions than in the cosmo- 
logical background and means that time variations in a 
will turn off in virialised overdensities even though they 
continue in the background universe. The fact that lo- 
cal a values 'freeze in' at virialisation, means we would 
observe no time or spatial variations in a on Earth, or 
elsewhere in our Galaxy, even though time-variations in 
a might still be occurring on extragalactic scales. For a 
cluster, the value of a today will be the value of a at 
the virialisation time of the cluster. We should observe 
significant differences in a only when comparing clusters 
which virialised at quite different redshifts. Differences 
will arise within the same bound system only if it has 
not reached viral equilibrium. Hence, variations in a us- 
ing geochemical methods could easily give a value that 
is 100 times smaller than is inferred from quasar spec- 
tra. We conclude that the consideration of the evolution 
of inhomogeneities, notably the one inside which we live, 
is essential if we are to make meaningful comparisons of 
different pieces of astronomical and terrestrial evidence 
for the constancy of a. 
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